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ABSTRACT 


The objective of this investigation is to derive the reliability and the 
associated probabilistic failure criterion for composite materials under 
combined stress. In the analytical derivation, the concept of joint probability 
was used and applied to the Weibull distribution function. In applications, given 
the experimental measurements of the necessary statistical parameters for the 
specific composite, the probabilistic criterion of the composite failure and the 
reliability of the specific structure can be predicted. 

Graphical representations for the joint reliability and joint failure contours 
were made in two and three dimensional space for the several different sets of 
Statistical strength parameters to illustrate the effect of parameters on reliability. 
such understanding will enhance selection of fiber and matrix (which have their 
own statistical strength parameters) and can lead to improvement in reliability of 
some composite components in an aircraft structure. These reliability and 
failure concepts can also be used in repair problems by selecting the proper 


composite material with the appropriate statistical parameters. 
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1. INTRODUCTION 


A. FAILURE AND RELIABILITY OF THE COMPOSITE MATERIAL 

In a composite, the matrix and reinforcement materials are combined to form 
a macroscopically multiphase material . The macroscopic state gives rise to a 
material which is internally redundant against local defects and damages. 
[Ref.1] The ply is the basic unit in constructing composites which consists of 
fibers embedded in a matrix. It is possible, therefore, to relate composite design 
properties to corresponding constituent properties. As a result, it is possible to 
design structural components made of fiber-composite when the design 
requirements and the properties of the candidate constituent materials are 
Known. For high performance structures, the unidirectional fiber-matrix 
composites consist of strong stiff fibers embedded in a comparatively ductile 
matrix . Fora rational structural design, the strength of such unidirectional ply 
needs to be expressed in terms of the stress tensor. Because the ply is 
composed of numerous fibers, the composite failure process is sequential, 
therefore ad hoc failure criteria can not be applied to analyze the composite 
failure. Furthermore, such conventional failure criteria such as Mises and 
Tresca's criteria are related only to a single parameter (mean function of 
failure). As a result, they can not be applied to multiple parameter functions, so 
Suitable general relations to represent failure contours of composite materials 
are required. 

In the application of a one-dimensional failure criterion, the magnitude of 
stress (one component) at each spatial location of the structure is mapped into 


respective points on the stress space (which is a line). If all the points within the 


domain are bounded by the tensile strength limit and the compressive strength 
limit, then the entire structure may be considered safe. The reliability under 
unidirectional stress, then, is the product of the reliability of each spatial 
segment at the respective stress level. Under planar combined stress, the 
magnitude of stresses (3 components) at each spatial location of the structure is 
mapped into respective points on the six dimensional stress hyperspace. So the 
stress that is interior to a domain bounded by the failure surface is safe and the 
failure criterion is the analytical, graphical, or numerical representation of the 
failure surface. To formulate the failure criterion, consistency of mathematical 
operations can be assured by adhering to the quotient rule of tensor variables. 
Historically, the formulation of many failure criteria were intuitively taken to be 
based on failure mechanisms but in fact, they are phenomenological because 
of the inconsistency between macro and microscopic failure modes. 
Phenomenological failure criterion may be considered as a mathematical 
model of the material transfer function relating the external excitation (stresses) 
to the material's response (failure). The phenomenological mathematical 
model is intended to aid experimental design to facilitate interpolation, 
correlation, and retrieval of experimental observations. They do not, in general, 
address the statistical variability of strength. 

For a combined stress failure mechanism, when one stress component does 
not affect the strength of another component, the failure is mechanistically 
uncoupled; conversely when one stress component affects the strength of 
another component, the failure is mechanistically coupled. In addition, if 
probabilistic failure is also taken into consideration, the condition for the failure 
mechanism is not fixed; it depends on the stochastic combinations of the 


intrinsic strengths. Under combined stress, mechanistic uncoupling with 


Statistic independence can interact to produce phenomenologically coupled 
contours. On the other hand, mechanistically coupling with statistical 
dependency can also produce phenomenologically coupled statistical contours. 
These cases can only be distinguished by examination of both the contour in 
the stress space and the statistical strength distribution in the probability space. 
In any phenomenological failure criterion, the shape of the failure envelope 
is never completely known until experiments are performed for all possible 
combined states of stress. As a result, complete experimental determination in 
two or more dimensions place an impractical demand on time and resources 
because of the large number of tests needed. However, with the establishment 
of a probabilistic failure criterion, this experimental failure contour can be 
related to the shape parameters (the variability) in the material’s principal 
direction. With such supplemental information, failure contours connecting the 
mean value can be established by a smaller number of samples. Therefore an 
analytical model can complement the failure characterization by reducing the 
required number of experimental measurements, providing bases for 
comparison of the material’s performance, and simplifying data and data base 
reduction. In structural analysis applications, the reliability of a specific structure 
under a combined stress can be determined by mapping the spatial domain into 
the stress domain through stress analysis. Using the probability failure criterion, 
the reliability of each spatial location corresponding to the respective stress 
tensor and magnitude can be calculated. The joint reliability of the individual 
Spatial location can in turn be combined to predict the reliability of the entire 


Structure. 


B. OBJECTIVE OVERVIEW 

The objective of this investigation is to derive an analytical model for 
composites failure and reliability under combined stress conditions. The 
Weibull distribution function was used in the probabilistic modeling. Chapter 2 
presents the theoretical formulation for failure and reliability, and the resulting 
equations are shown graphically in the following chapters. Chapter 3 shows 
the joint reliability and linearized failure Cumulative Distribution Function (CDF) 
in two and three dimensional space. The effects of parameters which depend 
on the material were analyzed. Chapter 4 shows the joint failure CDF in two 
and three dimensional space. The effects of material parameters were also 
analyzed. In Appendix B, these probability distributions were compared to 
experimental data availabe in the literature. 

Finally, Chapter 5 provides conclusions and some remarks for the future 


work in this area. 


Hl. TH AL FORM TION 


A. RANDOM VARIABLES 

A failure criterion for a solid prescribes the risk of failure in terms of the 
causes of the failure. Mechanical catastrophic rupture is considered herein. 
The effect of mechanical failure may be considered to be caused by either 
stress or stain depending on the microscopic failure processes. lf the 
constitutive relation of the material is available, conversion between stress and 
strain can be readily made and the two variables may be considered as 
equivalent. Similarly, a probabilistic failure criterion which describes the 
probable risk of failure can be expressed in terms of stress or strain as random 
variables. Stress is selected as the random variable herein with the 
understanding that the variable can be changed to strain with the appropriate 
constitutive relation. We denote the probabilistic random variable as Xj, for the 
scalar component of the stress tensor , and xj for the realized random variable 


(i.e. the strength X; under the respective stress component o}). 


B. PROBABILISTIC CONDITIONS OF FAILURE 

Under combined stress conditions, when more than one component of the 
stress tensor assumes non-zero values, the effect of each component of the 
stress tensor contributes to the failure of the composite. Failure of the 
composite does not occur if failure by each individual stress component does 
not occur. 

For the composite subjected to a planar state of stress, we denote the failure 


events as: 


A: Failure caused by normal stress along the fiber direction - oy 
B: Failure caused by normal stress perpendicular to the fiber direction - oo 
C: Failure caused by shear stress along the fiber direction - og. 
Then the failure of the composites is either by A or by B or by C or by all three as 
illustrated in Figure 2-1: 
P(A +B +C) = P(A) + P(B) + P(C) - {(P(AB) + P(BC) + P(CA)} + P(ABC) 
| (2-1) 
As a special case, if all three events are mutually exclusive, then Equation 
(2-1) can be simplified by 
P(A +B + C) = P(A) + P(B) + P(C) (2-2) 





Figure 2-1 : Joint Probability 


C. JOINT DISTRIBUTION FOR FAILURE 
1. Discrete Case 

In the development of the analytical expressions for the probability of 
failure, we first developed the general form of the distribution for three random 
variables using the notation Xj for the random variables and xj for the realized 
random variables. For the planar stress, the customary notation for the range of 
i = 1,2,6 is used. These general expressions are then specialized for the 
combined stress cases by replacing the realized random variables by oj. 

The probability of failure under the combined influence of three random 
variables can be expressed by the joint probability density function (pdf) of the 
individual realized random variables: 

FOS On Pe 04,X, = 6.) =1(0,,6,,06,) (2-3) 
where f(o1, Oo, dg) 2 0 and YD f(o1, 65, 6g) = 1. If random variable Xg is to 
take on any one of the values of 067, O62, 963:---, OG, then the joint pdf at 


some combination of X14=04;, X2=09); and Xg=0¢6; Can be expanded by 


P(X, = 6X) = On XQ = 6.) = F( FF 55,5 ¢y) (2-4) 
Nn 
P(A, = 04 Xq = oF) = (9, 53) oa 21 (F 4 Fo) F@y) (2-5) 
| 
P(X, = 55,X6 = Og) = 195; O64) = LA (O 455, © 6x) 
i=0 (2-6) 
mM 
P(A, = Fy Xe = G gy) =F (F454) = LF (O41 Fy» Fg) 
j=0 (2-7) 


and the marginal probability function can be obtained for the three random 


variable case: 


m on 
Co) aoe 22 (F927 F2 7 Pek) 


(2-8) 
yy 
PX = 6.) = 1 Oe eS {(0,.,0,.,0 
(X, 2) ( 2)) we > BF 6k) (2-9) 
1 m 
(Ag = Cap ie (oleae LL 4 (54:55; 64) 
i=0j=0 (2-10) 
which can be written 
ea n 
2F,(6,)) = 1, =f »(52;) = + Dif 6(F 64) = 
i=0 j=0 k=0 (2-11) 
[om 2A 
pap: LF (05:5 5/5 64) =a 
i= Oj= Ok=0 (2-12) 


2. Continuous Case 
lf three random variables are considered, then Equation (2-4) through 
(2-7) can be expanded by 
P(X, < One < Oni < 6.) = F(0,,0,,0,) 


me 2 ° 
= ee £ (uv,wdudvdw oan 


ee 
RCGsi =e) ncn J J _Jf(u,v,w)dwdudv - 


6 ° 2 as 
P(X, $0,,X_, 5 6,4) =F(o,,0,) = J J Jf (u,v,w)dudvdw P 


15) 
6 ° ie 
P(X, 50,,X, $ G,) =F(5,5,) = J J _Jf(u.v,w)dvdudw oi 
and the marginal distribution can be obtained by 
o oO oo 
P(X, < 6,) = F(o,) = J J JfWv,w)dvdwdu oa 


oe ema 
P(X, <6,) = F(o,) = ie J Jf(uv,w)dwdudv 


(2-18) 
Coe 
mes, < 6) = F(.) = J J _Jf(u.v,w)dudvdw ie 
which can be written 
f f {#(u,v,w)dudvdw= 1 
a (2-20) 


From the definition of conditional probability, the conditional probability of 


Oo given that 6; occurs can be represented by a function form: 
1(0,,0,) 

f(o,|0,) f4(9,) (2-21) 

where f(01, Oo) = P(X4=04, Xo=00) 
a. General Case 

If we consider the combined stress case, the random variables Xj, 
Xo, and Xg can be assumed to be the failure strength in the 1, 2, and 6 direction 
respectively. And under a combined stress condition, the failure occurs when 
either failure strength X71, Xo, or Xg is reached. Therefore the probability of 
failure under combined stress X1 and Xo can be represented by 

P(X1 + Xp) = P(X4) + P(Xa) - P(X1, Xo) 
and then the above relation can be related to the joint distribution function as 
follows 

F(04, O65) = P(X4 $04, X98 00) = P(X4+ X9; X14 $04, X20) (2-22) 
and similarly for three random variables case: 

F(61, 69, 6g—) = P(X1S01,Xos00,X6s06) =P(X1+X0+X_E | 
X1504,X2500,XGS0¢) (2-23) 


Equation (2-22) is proven in Appendix A-1 using the Weibull 
distribution function for the independent case. 

lf the joint distribution of failure under the combined stresses of 64, do, 
and og is considered, the strength failure criterion of a composite is the 
minimum of the failure strength X1, Xo, and Xg, that is, at Xp = X4AXoAXe. ( We 
are not dealing with the minimum magnitude among Xj, X9, and Xg but rather 
One value among Xj, Xo, and Xg which is less than the intrinsic failure 


strengthes of 61, 65, and 6g becomes the minimum value here.) So 
HO,05,0,6) = PIX, =X, aX, a X,OVET( 5,55, 4) } 


where ( 61, 09, 0g ) is the three dimensional stress domain. If the relations 
among 61, 69 and og are given by 

0; =B1209, 09 = Bogog, dg = Be10;4 (2-24) 
then, the joint pdf P(X1< 01,X2S00,Xg<0¢g) can be evaluated. If og is assumed 
to be zero, then the joint probability functions, P(X1<04,Xa<09) can be divided 
into two probability distribution functions corresponding to S; and So domains 
as shown tn Figure 2-2. The domain Sj is applied for the case when 
[X14 / Xo <Byo and Xj < 04] which the composite fails by failure strenth 
Xo = X44Xo, whereas the domain So is applied for the case when 
[X4/X2 2 By2 and X92 < og] which the composite fails by failure strength 
X41 = X4AXo9. So the following joint relation can be obtained: 

P(X $04, X2 < 69) = P(X1/Xo < B12 and Xj < 04) 

+ P(X4/Xo 2 B12 and X20 < 69) (2-25) 

Applying Equation (2-25) to domains described in Figure 2-2, then 

P(X, 0,,X, $05) = ah (u,v)dudv + J. | f(u,dudv 


fe 
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Using end limits and considering P(X4<o1,X9S09) = F(o4,09),then the joint 


distribution function for the general case can be represented by 


oe is 5 es 
F/o,,0,) = , t(u,v)dvdu + f (u, v)dudv 
(rte) 1, Laltowtentud, as (2-26) 


where u and v are dummy variables. 


NEI : 


LL 
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Figure 2-2 : Domains For The Joint Distribution 





b. Specific Case 
The joint probability of the bi-axial case (Equation (2-25)) can be 
expanded for specific cases. 
(1) Independent Case. If the events X; =6; and Xo =o are 


independent for all o; and oo, then the joint probability can be represented by 
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P(X4=074,X9=09) = P(X4=04) P(X2=09) 
or equivalently 

f(04,09) = f4(04) f2(09) (2-279 

If three independent random variables are considered: 

P(X4=04,X2=00,X6=06) = P(X1=04) P(Xa=09) P(Xg=06) 

f(54 69,6) = f(04) f(69) f(g) (2-28) 

Similarly, if the random variable X4 and X9 are independent for all 
Oy and oo: 

P(X4S04,X2S09) = P(X4S04) P(X2s09) 
or equivalently: 

F(04,09) = F4(64) Fo(o9) (2-29) 
where F4(o4) and Fo(o9) are the marginal distribution functions of X4 and Xo 
respectively. Similarly, for the three independent random variable case: 

F(X4<04,X2S090,X6S0¢8) = P(X4S04) P(XoS00) P(Xg<0¢6) 
or equivalently 

F(04,09,06) = F(o4) F(a) F(o¢) (2-30) 

So the joint distribution function can be expanded a step further for 
the independent case. 
From Equation (2-26) (refer to Appendix A-2) 
F(5,,55) = F.(95,) + F (55) 


o Oo 
1 U 2 
ms Fg \ou- J f (WF,(B,.v)dv 
0 12 0 (2-31a) 
Equation (2-31a) can be further expanded by integration by parts 


F(O,55) = F.(S,) on F (55) - AF (9,)F,(55) 


Oo 0 
1 u 
7S | FU f(g Bu - Bia) F(v)f,(B,.v)dv 


12 (2-31b) 
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Comparing Equations (2-31a) and (2-31b) reveals that the 
arguments of the pdf and CDF are exchanged so either Equation (2-31a) or 
(2-31b) can be used depending on simplicity of algebra. 

(a) Weibull Distribution Model. If the Weibull function is applied to 


the independent case, then Equation (2-31a) can be expanded by 


eonat-od (3) It oo) | 
tala ta t-te 
tet ah tLe 


By rearranging terms, the joint distribution function for the 








Weibull model can be obtained by 


veel (R LR) © -() “(eim) 


a 





Pa Vy be” uy? (Byou | 
“Jel eG) “GE Js 
o Bs \ Bp B. B, (2-32) 
(2) Independent And Identical Case. Equation (2-31a) can be further 


expanded if the identical case is considered. 





For the identical case, F(o) = F4(01) = Fo(o9) and equivalently f(o) 
= f4(01) =fo(05), so expanding Equation (2-31a) 


6 6 
1 2 
F(o,,0,) =F(o,) + Fo.) - FFB — Jou + f (V)F(B,.v)dv 
(752) =F(9,) +F(O,)- J FUG lout] TFB, 2¥) 
(2-33a) 
and similarly Equation (2-31b) can be represented by 
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F(04,09) = F(o1) + F(o9) - 2 F(o7) F(o9) 
| °, °° 
ae} pore. ner 


For the special case when B19 = 1, 6 = 61 = 65 and Equations 


(2-33b) 


(2-33a) and (2-33b) can be simplified to 


F( Sy O5) = Z1ACe)) 2 ft WF Udu 


and rearranging finally gives the joint distribution function for independent and 
identical case.(refer to Appendix A-3): 
F(64,05) = 2 F(o) - (F(a))@ (2-34) 
Equations (2-33a) and (2-33b) can be evaluated for specific 
models such as the Weibull distribution function. 
(a) Weibull Distribution Model. For the Weibull model, the pdf (f(x)) 


and cumulative density function (CDF: F(x)) can be represented by 


)= (818) “exo(- 3) (2-35) 
i alka exp|- (5) } (2-36) 


Substituting Equation (2-35) and (2-36) into Equation (2-33a) 


(Opp) = : o| - F)} : - on - (3) 
I, (SLR) f- (Bayt ex0[- (gy |b 
BBY o-6NIt-9- FPP 








By integrating and rearranging the terms, the joint distribution function of the 


Weibull model under bi-axial stress 6; and 65 can be simplified to (refer to 
Appendix A-4) 


Cor) es Us oF é F] as e87) 


D. MEAN FOR THE JOINT DISTRIBUTIONS 


a 


1. General Case 


If X; and Xo are two continuous random variables having the joint 
probability density function f(o1,09), the mean or expectation of X;and Xo can 


be obtained by: 


=a E(S,) = J J of(o,0,)do,00, (2-38a) 


H., = E(o 5) ~ J J 9, {(6,0,)do,d0, ee. 

Considering the stress, the bottom limits of the double integration are 
changed, so Equations (2-38a) and (2-39a) can be expanded specifically for 
the stress case which we are concerned with 


i 0 oa — 6.,0,\do0.do 
1 E(S,) of (F 452) ie ee (2-38b) 


M,= E(o 5) = J J of (57 %,)do,do, Bee 


Similarly if three random variables are considered, then Equations 


(2-38b) and (2-39b) can be expanded by 


oO  c& _ & 


i 0.f/6,,0,,0,)\do.do do 
1 J J J oto, 2 6) [2° 6 (2-40a) 


oo oc co 


oS SlCr Gr GAG dodo 
z J J J onto, 2 6) gna aes (2-40b) 


is 


H,= J J J 5 61(502154)05,d0,de, cual 
2. Specific Case 
a. Independent Case 
If the random variables X14 and Xo are independent for all o, and oo 
(04,09) = 4(04) fa(c9) 
and Equations (2-38b) and (2-39b) can be expanded by 


Hes lies G.\1./G) 06.66 
i=] J of,(9,)f(%2)40,00, oa 


Weal s ns 
2 J J 2" 1( 1) 2( 2) te (2-42a) 


taking the integral and noting that F4(e%) = Fo(e) = 1, then Equations (2-41a) 
and (2-42a) can be simplified to 


,= am o,\do 
1 J "4 1) 1 (2-41b) 


i= Omir o-\ao 
2= J Saf 2(52)d0, on 


Similarly, for the three random variable case: 


Co © _ & oo 


ee J. J J 5,1,(0,)f 2(0,)f ,(0)49,do,do, = J of,(9,)00, 
(2-43a) 

ee JJ J 9,,(0,)f(%2)! g(% 6 A990 06 = Jo,f 9(F2)99 5 
(2-43b) 

Le FJ J 066.) f2(%2) 3(Sg)40 do ,do, = Jog 5(F6) 9. 
(2-43c) 


Observing Equations (2-41b) through (2-43c), it can be concluded 


is 


that the mean values of the random variables X;, Xo, and Xg are dependent 
only on their corresponding probability density functions. 
b. Independent And Identical Case 

For the independent and identical case, the joint probability function 
can be obtained by 

f(04,69,0g) = (04) f(o5) f(G¢) 
where f(o) = fy (0) = fo(o) = fg(o) 

Applying these relations to Equations (2-41b) through (2-43c), the 
mean values for the joint distributions can be simplified for the independent and 


identical case. 


Wi =p,=u = | of (o)do 
1 2 ay. (2-44) 


(1) Weibull Distribution Model. If the Weibull probability function is 


considered, Equation (2-44) can be expanded as follows 
iZ ayo” o\" 
= == | | eas SIS) alo) 
H= J (BB) | (5) | 
letting = (o/B)4, then 
co 
u=] By" exp( - y)dy 


Using the definition of the Gamma function, we finally get the 


mean function: 


aa prt ss = (2-45) 


E. JOINT RELIABILITY 
ile n f Reliabili 


The common notion of reliability is the confidence in the ability of a 


1 


device to perform adequately over a range of specified conditions. Such 
qualitative measurement of confidence and adequacy can be quantitatively 
calculated by probability. !n the formalism of probability, the conditions over 
which the device is intended to perform is defined as the random variable X: 
and the probability over the range of conditions is defined as the CDF over the 
range that the random variable experiences in service. That is, over the sample 
space o experienced in service, the probability of failure is defined by the CDF 
over the sample space: 
PCs Gil hie, 


Since reliability is the compliment of failure, it is defined by 


miG) = 1 FG) (2-46) 
2. Joint Reliability 


lf X4 and Xo are two random variable, then the joint reliability of X4 and 
X92 can be defined by 
R(o4,02) = 1 - P(X4<071,X9<09) (2-47) 
or equivalently for the three random variable case 
R(o4,02,06) = 1 - P(X1<04,Xo<09,X6<06) 
For the discrete case, equation (2-47) can be expanded by: 


R905) =1- >» df (uv) 


USO {¥25>5 (2-48) 


R65, og =1-e S'f(u,v,w) 


USO ;Vs9 one SO, 


lf the continuous case is considered, the joint reliability function of X; and 


X9 Is defined by 


R6,5,) = =1- i" yi (u,v)dvdu (2-49) 


and from Equation (2-49) the marginal reliability function can be obtained: 
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R(9,) =1- J _ Jf (uvydvdu pepe 


oo 
R.(o,)=1- f(u,v)dudv 
2(%2) J Jr (2-50b) 
lf the three random variable cases are considered, Equation (2-50a) and 


(2-50b) can be further expanded: 


0 oe) foe) 
Ryo, =1-f J Jf uv.wydvdwdu 


(2-51a) 

R,(5,) =1- J J J (u,v, wdwdud i 
OPEN ies 

R,(5,4) =1- J J Jf(tv,wdudvaw — 


R7(07), Ro(o2), and Rg(o¢) are also equivalent to the one dimensional 


reliability function which are defined by 


R(9,) =1- J f,wu 


(2-52a) 

F(a) els J Aen (2-52b) 
ar 

Re(Fe) = 1- J egiee (2-52c) 


a. General Case 
The general equations for reliability can now be applied to the 
calculation of the mechanical strength reliability of composites under combined 
stresses by specifying the random variables as X; and the realized random 
variables as oj. The joint reliability under combined stress o; and oo can be 


obtained by 
(6,55) =1- F(5,,05) = |- PIX=x, A X ,OVEr(O,,5,) } 
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where Xj and Xo are the failure strength in each coordinate(see Figure 2-2). 
Noting that the joint pdf under combined stress oj and oo is 
represented by the double integration of f(o1, o9) over the domains Sj and So 
and that the double integration of f(o], o2) over the domain Sg (see Figure 2-2) 
is defined as the joint reliability under combined stress 6, and oo, then the 


following relation can be obtained: 


636,\=1—-F(o,c = : anncuck: 
Ara, = 1-F(oyo)= J, J, oa 


Comparing Equation (2-53) to Equation (2-26), Equation (2-53) can 
be represented in another form: 


5,65) =1- [ Je tuvpdvdus | fiw, vidum 
12 ary (2-54) 
If the joint reliability under three combined stress is considered, 


equation (2-53) can be expanded by 


6.0.60. \= 1 —F/6, 6 =o.) = ‘i “fv wielewday 
AO Fy6)=1- Foy) =), Jo Ie oa 


b. Specific Case 
The joint reliability of bi-axial stress case can be further expanded for 
the specific case as in the probability distribution function. 
(1) Independent Case. If the random variable X; and Xo are 


independent for all oy and o9, then Equation (2-53) can be expanded: 


F(O,55) = i: JA v)dudv 


(2-56) 
After integration, Equation (2-56) can be simplified to (refer to 
Appendix A-5) 
R(o4,02) = {1-F4(04)}{1-Fa(oa)} (2-57) 
noting that Ry (o4) = 1-F4(o4) and Ro(o9) = 1-F9(o9) 
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R(o4,062) = Ry(o4) Ra(o2) (2-58) 
For the three random variable case, Equations (2-56), (2-57), and 


(2-58) can be obtained by 


R( 055,04) = a j. jt, f,(v) f, (w)dudvdw 


= (1-F4) (1-F2) (1-F6) 
= R4(01) Ro(o2) Re(o¢) (2-59) 

In application to a composite, if the probability of failure model (Fj) 
is Known, then the reliability of the composite can be calculated. The probability 
of failure model can be inferred from the physical consideration of the failure 
processes. 

When a filament composite is loaded along the fiber direction, 
local failure begins when the weakest fiber fails. The load carried by the broken 
fiber is transferred to the neighboring fibers. Upon additional increase in load, 
additional fibers fail, leading to the increase of failure sites distributed over the 
composite. The higher the load, the higher the density of such failure sites and 
the higher the probability of clustering. The spatial clustering of the fiber failure 
sites leads to stress concentration and ultimately causes the catastrophic failure 
of the composite. Harlow and Phoenix investigated the probabilistic modeling 
of the above sequential failure events and arrived at a modified weakest link 
model in which the link of the chain is a bundle. [Ref. 2] For a limited range of 
the random variable, it can be approximated by the two parameter Weibull 


model [Ref. 3]: 


noose -(8) | 


When a filament composite is loaded perpendicular to the fiber 


ZA 


direction, the failure mode is similar to that described by classical fracture 
mechanics. The exception is that for a composite, all inherent flaws are not 
randomly oriented, but are aligned with and propagate along the fiber. In such 
a case, the largest crack dominates and it is effectively a weakest link process 
over the physical volume. On these physical grounds we may use the Weibull 


model for both the transverse strength and shear strength: 


neoeont (8) 
reore (i) | 


Equation (2-57) can be specified for these Weibull distribution 


e122) -209-(2) -($2) | 
B, B (2-60) 


equivalently for the three random variable case: 


e090) =| -() af yb | 
Z Pe 6 (2-61) 


If we note F(o71, 69) = 1 - R(o4, 69), then the joint distribution 


functions: 


function obtained in Equation (2-32) can be simplified by 


- 1 ol i) 2) | 

GG.) = le ea een | ee 

(i) B, B, (2-62) 
-+-exp}-(%1) (22) (28) 

A (rplay sy) = 29 a i 5 | (2-63) 


(2) Independent And Identical e. If the random variables Xj and 


X92 are independent for all o, and oo and the probability distribution function for 
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each random variable is the same, then F(o) = Fj(o7) = Fo(o9). Similarly R(o) = 
R(o4) = Ro(o9) and the Equation (2-58) can be simplified to 

R(o4, 62) = R(o4) R(o2) (2-64) 

As a specific case, if the Weibull distribution function is 


considered, Equation (2-64) can be further expanded: 


RO 55) = on : @) z () | (2-65) 


For Byo =1, then o = 01 =09 and the joint reliability function for 


Weibull model is obtained. 


ini exp|- 2g) } (2-66) 


For the three random variable case, Equations (2-64) through 
(2-66) can be expanded by 
R(04,02,0g) = R(o1) R(o2) R(o¢) (2-67) 


(G5 5,5 5) = o10|- (2) 7 (52) 7 (52) | (2-68) 


R(o,0,60) = exp|- (3) | (2-69) 
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III. METRIC REPRESENTATION OF JOINT RELIABILITY UNDE 
COMBINED STRESS 


A. BACKGROUND FOR GEOMETRIC REPRESENTATION 

CDF and pdf are useful functions which contain all the relevant information 
about the statistical properties of a random variable, in our case, the strength of 
composite. These functions are required in reliability analysis, design, 
acceptance, maintenance, and operational logistics. 

In applications, in order to map out the statistical failure surface even for the 
bivariate case (the biaxial combined stress), an exceedingly large number of 
experiments are required. A large number of experiments are frequently 
impractical due to economic and time constrains and other considerations. 
Therefore, it is important to visualize the shape of the failure surface in order to 
narrow the range of experiments to be focused on the critical regions in the 
stress domain. 

In this chapter we present the geometric representations for the reliability 
functions for the independent case derived in Chapter Il. The graphical 
presentation is based on the probability plots and the failure surface 
representations, which were also investigated. Examination of these graphical 
representations will shed light on the appropnate experiments necessary for the 
identification of whether the failure processes are independent or dependent. 
When the independence is established, the entire failure surface (all 
permutations and combinations of combined stress case) can be calculated 


from the uniaxial strength statistics. 
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B. JOINT RELIABILITY FUNCTION 


1. Two Ran Variabl 
From Equations (2-56) and (2-66), the reliability function for the Weibull 


model can be obtained by 


“eg exp}: (8) | (3-1) 


RELIABILITY 


Ee alpiic=3 
@ aloha=25 


R(sigma) 





sigma/beta 


Figure 3-1 : Effects Of a In Reliability 


ZS 


If we consider o/B as a normalized variable, then the reliability function 
versus o/B can be plotted for various values of a. Figure 3-1 shows two different 
graphes for large and small values of a. For the small values of a, the graph 
shows a relatively smooth curve, but as @ increases, the graph approaches to 
the step function and these relations can be further expanded for the joint 
reliability problem. If the independent combined stress cases are considered, 
the joint reliability of the composite under combined stress can be obtained from 
Equation (2-58). 

R(o4,02) = R4(o4) Ra(oa) 
Expanding R(oj, 09) 
) Le) 
(ra) ENP \ | B, (2-60) 


where 


a 104 
1 Ps 


aa 25 
R(o,) =e x1) | R(O,) =exp|- Fe 
lf we assume that 
By = V12 B2 (3-2) 
and noting that 
0, =Byi2 02 
then the joint reliability for the independent case can be represented by 


R(o4,02) = R4(01) Ro(o2) 
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ot -(5) jot (aa 
; B, Bib, (3-3) 


If we further assume Cj9=V19/B 19, then Equation (3-3) can be simplified 
by 
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0102) -014-(2) Jof-fof2} | os 


From Equation (3-4), the reliability functions R7(o71) and Ro(o2) can be 
plotted independently in two dimensional space provided that a1, ao, and C190 
are given. Then Ry(o7) and Ro(o2) are combined to yield the joint reliability 
function of R(o1,09) as shown Figure 3-2. Observing the Figure 3-2, it can be 
noted that the joint reliability is affected by the smaller value between R4(04) 
and Ro(o9) for a given value of 04/84, that is, the value of the joint reliability is 
always close to the smaller value between R4j(o7) and Ro(o9). When the values 
of both R1(o71) and Ro(oo) are close to '1’ for the specific range of o7/B1, the 
value of the joint reliability is also close to 1, but as either Rj(o7) or Ro(o9) 
decreases, the joint reliability also decrease depending on the smaller value 
between Rj and Ro. We can also observe the effect of Cz for the given a1 and 
Q9. 

As shown in Figures 3-2a and 3-2c, the Ro(o9) shifts to the right from the 
R71(01) as C19 decreases from ’1’ whereas Ro(o9) shifts to the left as Cz9 
increases from ’1’ while Rj(o7) remains constant. So, the joint reliability follows 
the smaller function between Ry(o7) and Ro(o2) and appears to almost overlap 
when Cjo9 is far from ‘1’. As the value of C19 approaches to ’1’, the values of 
R4(o4) and Ro(o2) approach each other and when Cj9=1, R4(0;) and Ro(o9) 
cross each other at o4/B,=1. The joint reliability function exists to the left of 
R7(o4) and Ro(o9) and has a weaker reliability than R1(o7) or Ro(o2). So the 


joint reliability is always less than other uni-axial directional reliabilities. 
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Figure 3-2a : Reliability Vs 0)/B, (C12=0.5) 
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Figure 3-2b : Reliability Vs o;/B, (C42=1) 
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Figure 3-2c : Reliability Vs o)/B; (C12=5) 
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Here we introduce the concept of linearized failure CDF to show the 
effect of the joint failure CDF at the tail area. 

To get the linearized failure CDF, denoted by F*(o), use logarithmic 
algebra. If we take the logarithm of Rj(o7) and note that Ry(o7) = 1-F7(07), 


then 


a 
1 


Multiplying both side by (-1) and taking the logarithm again, the 


(3-5) 


linearized failure CDF can then be defined as 
* a a 
BF Gy) = In InfR(o,) }| (3-6) 
Substituting Equation (3-5) into Equation (3-6), then 
sf a 
Po (S,) = a In B, 


By the same way, linearized failure CDF corresponding to oo can be 


(3-7) 


obtained: 
sk °1 
Bs (S,) = a, In(C ») +In B 


Comparing Equation (3-8) to (3-7), it can be noted that Fo*(o9) is 


(3-8) 


vertically shifted and rotated when compared to Fy"(oj). This result is 
represented in Figure 3-4. 


similarly F"(o4,069) can be obtained by 


a Oo 


] ys 
F*(5,,5,) =In | + cg 
P (3-9) 


Comparing Equation (3-9) to (3-7) and (3-8), it can be noted that 


F*(04,09) is no longer a linear relation in terms of In(o7/B1). Figures 3-4a, 3-45, 
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and 3-4c show the linearized failure CDF for different values of C12. As shown 
in these graphs, Fo*(o9) is shifted upward and rotated clockwise as O12 
increases and this causes a higher probability of failure. The joint linearized 
failure CDF, F*(oj,09), then follows the higher value between F4"(o4) and 
Fo*(oo) and does not show the linear line, especially near the area where 
F4*(o4) and Fo*(o9) intersect each other. So we can observe that , as the value 
of C19 increases for the specified a; and a9, Fo*(o2) shifts upward and rotates 
clockwise causing the intersection between F y*(o1) and Fa*(o2) to go 


downward, or in other words creating a higher probability of failure. 
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Figure 3-4a : Linearized Failure CDF (C4 9=0.5) 
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Figure 3-4b : Linearized Failure CDF (Cj9=1) 
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Figure 3-4c : Linearized Failure CDF (Cj9=5) 
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As a special case, if a=a1=a9, then the joint reliability function and 


Equations (3-7) through (3-9) can be simplified by 


R(o,,6,) = cx a (1 + “)| ofl 
oF 
mr) = aa (3-11) 
(om) = afin(C,,) +] eh 
1 (3-12) 
Sod Co alate | +In(1+C_,) _ 
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Figure 3-5a : Reliability Vs 0/3) (a@,=a2=5) 
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Figure 3-56 : Reliability Vs o4/B, (a@=a=15) 
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Comparing Equations (3-12) and (3-13) to Equation (3-10), it can be 
noted that Fo*(o9) is shifted vertically by a In(vj90/B19). F*(o71,09) is also shifted 
vertically by In{1+(v40/B;92)} regardless of values of C19. Figures 3-5 though 
3-7 show the reliability in two and three dimensional space and linearized 


failure CDF in two dimensional space for a = a1 = Go. 
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Figure 3-6a : Joint Reliability In 3-D (a,=a9=5) 
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Figure 3-6b : Joint Reliability In 3-D (a;=a 
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As shown in Figure 3-7, each failure CDF is parallel to each other and 
the joint linearized failure CDF almost coincides with Fy*(o4) or Fo*(o9), 
depending on the value of C19 and is parallel to both F4*(o4) and Fo*(o9). So 
in the physical sense, C15 iS very important for estimating the failure and 
reliability. The joint reliability is affected by the weakest reliability function and 


this is true for the joint linearized failure CDF. 
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Figure 3-7a : Linearized Failure CDF (a;=a9=5) 
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Figure 3-7b: Linearized Failure CDF (a,=a7=15) 
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2. [Three Random Variable Case 
lf the three random variable case is considered, the joint reliability under 
the combined stress can be obtained from the Equation (2-79): 
R(04, Oo, 0g) = Ry(07) Ro(oo) Re(o¢) 
Expanding this equation for the Weibull model: 


ro14284)-29|-(2) -(52) -(22)] 
1 2 6 (2-61) 
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Figure 3-8a : Reliability Vs o)/B, (C12=0.5, C16=1) 
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Figure 3-8b : Reliability Vs o)/B, (C12=C16=1) 


43 


lf we assume that 
By = V12 Bo, Bo = V16 Be 
and noting that 
61 = By9 O92, 1 = Big 5g 


where 


then the joint reliability for the independent case can be represented by 


0.29 


Figure 3-9 : Joint Reliability In 3-D (C16=1) 
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(3-14) 


(3-15) 


(3-16) 


R(61,67,06)=Ry (01)Ro(o7)R6E(o6) 


Oo 


of fC 
] 


O 


) Jol-[et) Jor le) 
=exp}-|——| fexp}—-|C.=-| fexp}-|C = 
B, oa Ne (3-17) 


C sralkenc sss 
eae = en ae = (3-18) 


where 


Furthermore, if Cyo and Cy6 are given, then Cog can be calculated 


using the chain rule: 


V B 


econ 
Cog RB 


©) 


0, Bp B, peg 16 


59, 8 B. 


2 Be jo 2B ig 12 (3-19) 
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Se 
oO) 
| a 


_ 
©) 


SO, if 1,9, a6, Cy2, and C16 are given, then the joint reliability 
function, R(o1,09,0g), can be plotted in two dimensional space with respect to 
01/4. Here we can also introduce the concept of linearized failure CDF to show 
the effect of the failure CDF at the tail area. To get the linearized failure CDF, 
denoted by F*(o), use the same procedure described in Equations (3-5) through 


(3-9), then 


Ra (°,) =e! 5] 
| (3-7) 


1 


Ra (55) =aIM(C,,) +| ial 


1 


Ba (ae) =a) tS alia (a 


(3-8) 


1 (3-20) 
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a a 
2 


F* -in(2 ia (ces) | 
ee J fF: C1252) Coy (3-21) 


So it can be noted that the joint linearized failure CDF, F"(o1,09,0¢) does 
not show a linear relation in terms of In(o4/By) except when a4=a9=G¢, instead 
it is affected by F4"(o1), Fo"(oo), and Fg*(og). That is to say, F*(o4,09,0¢) is 
dominated by the weakest value among F1*(61), Fo"(o9), and Fg*(c¢) so 
F*(64,09,0¢) is always Shown on the top of each of graphs as shown in Figures 
3-10a and 3-10b. Figures 3-8a through 3-8b show the reliability for different 
values of C19 when aj, ao, C16 is fixed. These figures are expanded to three 


dimensional space as shown in Figure 3-9. 
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Figure 3-10a : Linearized Failure CDF (C19=0.5, C1g=1) 
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Figure 3-10b : Linearized Failure CDF (Cj9=C46=1) 


47 


As a special case, when a=a1=07=G¢, then the joint reliability and 


linearized failure CDF functions can be simplified by 


F(Fp5 9,5 5) = ox et (1+C,,+ c*)| 
(3-22) 


0} 
F.(6,) = of1N(6,.) ' ora = 


re) 
A = a a 4 
F 5 ,,05:0,) =alr{ 5t]+In(1+ C4, +54) 


Observing Equations (3-22) through (3-24), it can be noted that Fyj*(o7), 


(3-24) 


Fo"(o9), Fe*(og6), and F*(o4,02,0¢6) are parallel each other. And the Joint 
linearized failure CDF, F*(o1,00,06), is also affected by the largest function 
among Fy", Fo*, and Fg*. The parameter a is related to the rotation whereas 


C12 and C16 are related to the vertical shift of the graph. 


48 


ROBABILISTIC FAILUR il NDER MBINED STR 


A. PROBABILISTIC FAILURE 

For a combined stress failure mechanism, when one stress component does 
not affect the strength of another component, failure is mechanistically 
uncoupled, as shown in Figure 4-1. When one stress component affects the 
strength of another component, failure is then mechanistically coupled, as 


shown in Figure 4-2. 
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Figure 4-1 : Mechanistically Uncoupled Mechanism 
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Figure 4-2 : Mechanistically Coupled Mechanism 


lf the probablistic failure is considered, the mechanistically uncoupled case 
may contribute to statistical coupling so both failure mechanisms are not 
interacting and the probablistic failure mechanism is distinguished from the 
deterministic case. In order words, the condition for the failure mechanism is 
not fixed, it depends on the stochastic combinations of the intrinsic strengths. 
So, for the combined stress case, uncoupled independent statistical effects 
produce phenomenologically coupled statistical contours and the mechanically 
coupled dependent statistical effects also produce phenomenologically coupled 


Statistical contours. In the following subsection, the joint probabilistic failure 


90 


CDF will be discussed for the two and three random variable cases and the 


effects of the parameters will be shown graphically. 


B. Two Random Variable Case 
To show the joint failure probability, we here introduce the Weibull 
distribution function. From Equation (2-62), the joint failure distribution function 


can be obtained for the two random variable case. 


a QO 
2 


F(O, 5) =1- a | ~ ts 
B, Bs (2-62) 


From the above equation, (64/84) and (69/B5) can be computed if the values 
of F(o4,05) and C;9 are given. Figure 4-3 shows the joint failure CDF in three 
dimensional space and Figure 4-4 shows corresponding joint failure contour in 
two dimensional space for the specific value of F. As shown in Figures 4-3 and 
4-4, the shapes of the failure function and failure contour depend on the 
parameter a, which is material dependent. For a small value of a, the failure 
contour shows a smooth curve but as the value of @ increases, the joint failure 
contour approaches to the shape of a rectangle. This phenomenon is the same 
for the joint reliability contour as shown in Figure 4-5. It can also be noted that 
the failure contour is affected by a large value of a between a, anda». So if 
anyone of the values of a@ is larger when compared to another value of a, then 
the failure contour approaches to the shape of a rectangle, depending on the 
magnitude of a. These relations are shown in Figures 4-6 and 4-7. Another 
result, shown in these failure contour graphs, is that the distance between the 
contour line is small when a is large compared to those when ais small. So the 


intersection of the function on both axis depends on the value of each a. 
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Figure 4-3a : Joint Failure CDF In 3-D (a;=a9=5) 


D2 


Joint Fartlure CDF 


0.87 


0.58 


0.29 


0.00 


XY) 
ee 


y 
\) 
re 
WKN 
WS 
WW’ 
WwW 
AK 
Mi Me 
OXK 
X 
Wie 
Weg 
RAHI 
WKY 
I XKNXKRK 
Nes 
AX 
A 
A 
WK 
i 
Wn 
WX 
MAK 
MN 
MA 
AX 
ANN 
WW 
© 
{ 


ny 

\) 

\ 
A 
‘4 
\ 
Wy 
NN 
WW 
WN 
Mi) 
WN 
Wy 
\ 
vc 
\ 

| 


‘" 
4, 
( 
\ 
s 
bOI 
\\ 
My 
AVN 
MMA 
My 
NX 
WV 
MA 
in 
WN 
VM 
WA 
Nv 
M\\ 
WV 
oe 
0 XX 
ty 
9 
L 
§ 
y, 
; 


4 


\ 
Ys 
\ 


; 
é 


Q 
§ 
0 


Figure 4-3b : Joint Failure CDF In 3-D (a@)=a=25) 
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Figure 4-4a : Joint Failure Contour (a~1=0=5) 
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Figure 4-4b : Joint Failure Contour (a)=09=25) 
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Figure 4-5a : Joint Reliability Contour (a;=a5=5) 


36 


1.00 


0.50 





0.00 0.50 1.00 


Figure 4-5b : Joint Reliability Contour («;=a.9=25) 
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The contour line of F(a4,09) =0.5 exists around 1.0 in both the normalized x 
and y axis regardless of a, but as a increases, the F=0.5 contour line 


approaches to 1.0 in both normalized axis. 
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Figure 4-6a : Joint Failure CDF (a,}=5,a9=3) 
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Figure 4-6b : Joint Failure CDF (a )=25,a9=3) 
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As an application, if we compute the stress of a specific part of system such 
as a rudder or ailerons using the finite element method, then we can estimate 
the probabilistic failure and reliability of that specific part of the system. 

The parameters a and B, which were used in the failure CDF and reliability 
function, are material dependent so these values should be determined through 
experimentation and then these parameters can be applied for that specific 


material. 
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Figure 4-7a : Joint Failure Contour (a@,=5,a7=3) 
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Figure 4-7b : Joint Failure Contour (a =25,05=3) 
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C. THREE RANDOM VARIABLE CASE 


If the three random variable case is considered, the joint failure CDF can be 


developed using Equation (2-63). 
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Figure 4-8a : Joint Failure CDF (C46=1) 
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Figure 4-8b : Joint Failure CDF (C;¢=3) 
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From the above equation, normalized axis coordinates can be computed if 
the values of F(o1, 09, dg), C12, and C16 are given. And if we specify the value 
of C76, then we can plot the joint failure CDF in three dimensional space with 
respect to 04/4 and 65/4 as shown in Figure 4-8. As the value of C76 
increases, the possibility of failure increases, that is to say, the intersection on 
the normalized x-axis moves to the left because d¢/Bg term is embedded in the 
64/B4 axis. Figure 4-9 shows the corresponding failure contour line in two 


dimensional space for the specific value of F(o4, 69, dg). 
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Figure 4-9a : Joint Failure Contour (C16=1) 
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As proven graphically in the three random variable case, the probability of 
failure increases as the number of random variables increase. So in a physical 
sense, a specific composite material will show the highest reliability when the 
external loads are applied in the pricipal axis direction because there are only 
two random variable, X; and Xo, in the reliability fuction. But when the external 
loads are applied in the off-principal axis direction, there exist some shear force 


random variable, Xg, which decreases the reliability of the material. 
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Fig. 4-9b : Joint Failure Contour (Cj90=3) 
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As an application, if the stress of the maneuvering aircraft is computed in an 
arbitrary direction, then the stress can be transformed into fiber, shear, and 
transverse matrix directions which can then be used to predict the probability of 
failure at that specific load. Sample calculations were made and the results 


were analyzed with respect to experimental data in Appendix B. 
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V. CONCLUSIONS 


This study was directed towards deriving the joint reliability and joint failure 
CDF under a combined stress condition. The results were specialized for the 
Weibull distribution function based on the observations that the composite 
material failure is adequately represented by the weakest link model. The effect 
of the statistical strength parameters, a and B, which are material dependent are 
illustrated using two and three dimensional graphical representations. 

To analyze the reliability of the composite, the inherent statistical strength 
parameters in the composite's principal direction need to be experimentally 
measured. Substitution of these parameters in the probabilistic failure criterion 
will allow for the estimation of the reliability of the composite for any state of 
stress. 

Comparison of joint failure distribution under the restriction of independence 
to experimental data suggests that mechanistic coupling of the failure 
mechanism needs to be included in future extensions of the formulation and 
data of much larger number of samples has to be performed at critical stress 
ratios to conclusively examine probabilistic independence. 

Further studies may be extended to apply these reliability and failure CDF to 
the specific part of a system which is made of a composite material. To do this, 
the stress of the specific part of a system such as an elevator or ailerons should 
be analyzed using such as finite element method. These equations can also be 
used in a repair problem which requires the least cost and most effective 


method to analyze the combined stress in a specific part of a system. 
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So, if we analyze the combined stress, the appropriate fiber and matrix with 


proper parameters can be selected to make the composite fit to that specific 


part. 


68 


APPENDIX A 
Chapter II Equations 


1. For the independent case: 
P(X1+Xo : X4S<o04, XoS09) = P(X1<01) + P(Xos09) - P(X4<01) P(Xo<09) 
= F4(64) + Fa(oa) - F4(04) Fa(oa) 
If the Weibull distribution function is considered, F4(o4) and Fo(o9) can be 


substituted into the Weibull distribution function: 


esrensemsseref-en-(8) Ibo (8) | 
-b-est-(31) fh-ot-(54) | 


By expanding and simplifying 


a a 
0; | ge 
i nieS = O75 = 0) = 1— XP 65 aa 


which is the same result as the Equation (2-32) 


i 


2. For the independent case: 
F(O,0,) = ie [fn 


Taking the integrals inside the parenthesis first and substituting limits upon 


a (u)f A yay 


f(v imanlaus J f { i 


12 


evaluation: 


F(0,,0,) = [ofr ec) — FB. |e 
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rT f (V)LF (ce) - F(B,,v) }dv 
=f f (u)du, F, ( )=f t lu 


Noting Fy (ce) = Fo(ce) = 1 from the definition of CDF and expanding 


F(o,0,) =f f (udu - jf (OF {3 


a 
“f4 svidy — J f .(v)F,(B, ,v)dv 
Taking integrals and noting that F7(0) = Fo(0) =0 


F(0,0,) = F(S,) ~ F (Fo) - J FFB |e 


le 


16} 
2 
“s TMF (BV )av 
3. From Equation (2-33a) 


Foy 55) = 2F(o) - 2 ft Fd 


By integration by parts 
Foy 55) = 2F(o) - 2HF(o)) - JFuFtydu 


noting 
i (u) F(u)du = 0.5{F(o)}* 


then the joint distribution function for independent case can be obtained by 


F(o4,09) = 2 F(0) - {F(a)}? 
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4. By expanding and integrating 


(ore) [1-e0e|-(5 | +09 -F) | 











By substituting and expanding 
(2492) 0|-(5) |] +1 enol-(8) 
5,4) B 5, a S, a 
+|oxo|-(2) [-a}-[ se Jol - , -[as5| t-1 
Ce} etalenl-)-Ce9)] 
+o cle (stow i) Io 
By rearranging terms and simplifying 
B’ o,\ Gay 
me. 12 ee ae 1 
eves) Basia = (a5) 
(gst pa re B 


noting 61 = By209, then the joint distribution function of the Weibull model 























under bi-axial stress 6, and oo can be simplified by 


fi] 


a 
10) 


F(5,,55) =e ox0)- (=) 


a 
(8) 
B 
5. From Equation (2-50) 


(0,05) = J. J fut s(v)dudv 


E J, FA oo) - Fo.) Hf (v)dv 
noting Fy (ce) = Fo(cc) = 1 
Fo, 05) = 4 - F (0) 3] f ,(v)dv 


integrating then we finally got the Equation (2-57): 
R(o4,69) = {1-F4(o)}H1-F2(o)} 
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APPENDIX B 


Structural Applications 


In many structural applications for aircraft, the composite laminae are 
oriented at an angle to the tensile load. Take for example the skin on a wing, 
where the filament angle is oriented primarily along the direction of tension due 
to bending. Transforming the stress by an angle @ to an off-axis introduces a 
shear coupling that will compensate for the bending/rotational coupling of a 
swept wing. Such a material loading relation is representative of that indicated 
in Figure B-1. We will explore the combined stress state for such a 
configuration and the structural reliability as affected by the off-axis angle 0. 
Since we are addressing a Spatially homogeneous state of stress, both stress 
analysis and the substitution of the combined stress in the probabilistic failure 


criterion can be carried out explicitly. 


Xo 





Figure B-1 : Transformation Of Coordinates 


is 


For the stress analysis, a load p is applied to the composite material and we 
denote the direction of p to be x1', then p can be transformed to the principal 
direction of the composite; x1 and xo where 

X1 : fiber direction 

x9 : transverse matrix direction 
and the angle between x1' and x1-axis can be defined as 8 which is positive 
clockwise. Using the tensor relation, the stress can be transformed from x4'-x9' 
axis to X4-x9 axis as follows: 

Oij = Ami Anj Omn’ (B-1) 
where 

l,j, m,n = 1, 2 for 2-dimension 

oF a 

77 ON) [Opie 

aij = /COS(x1’, Xe COS, ind = ae -sin(®) 


cos(x9', X14) cos(x9',x9) sin(®) cos(@) (B-2) 
Expanding Equation (B-1) for 0) ;, 6)9, 699, then 
011 =4m1 4n1 Omn 
= 11 411 911 +211 221 9j9 +42) 81205) 7 apa 
= a21*o11' 
01] =p cos2(9) (B-3) 
For the same reason 
012 =-p cos(8) sin(@) (B—4) 
077 =P sin(0) (B-5) 


To calculate the reliability associated with the applied load p, the stress 
components in terms of p (Equations (B-3), (B-4), and (B-5)) can be substituted 


into Equation (2-83) to find the joint failure CDF for the specified external load p. 
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We note that Equation (2-83) is specialized for the case that the effect of oj, 


617, and 079 On the respective failure modes are independent. Substitution 


yields: 


a a 
a 
1 2 G 


F(p,6) = 1- exp} — peos*(6)) _/(psin*(e)) _ (pcos()sin(@) 
B, B,, B 


where 
alte oll) °2 = 222726 — 912 
Comparing Equation (B-6) to (2-83) and (3-18), it can be noted that 
B10 =cot@(0), B16 = -cot(6) 
C12 = V12 tan@(@), C16 = -V1¢ tan(6) 
This joint failure CDF in terms of p and 6 can be plotted in terms of p for an 
specified 6. Similarly, from the Equation (B-6), the linearized failure CDF can be 


obtained by 


0 § a 
1 2 o. . 


2 ae . 
7 pcos‘ (6) psin*« (6) pcos (6)sin (6) 
cca eos cd Ua ce 


(B-7) 
In order to investigate the dependency of the failure mechanism, 
experimental data from Sun and Yamada [Ref. 4] is examined. From the 
experiment, the following parameters were obtained for fiber, transverse, and 
shear force using uni-axial tests. 
a, = 20.5, B) = 127000 
a7 = 5.66, By = 1070 
ag = 8.96, Bg = 3078 
Substituting these strength parameters into Equations (B-6) and (B-7), the 


le 


probability of failure for load p applied to different angles can be examined. 
Figures (B-2a) through (B-6b) compare the joint failure CDF to the experimental 
data for different angles. These data will be examined individually for each of 


the five angles where experimental data is available. 
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Figure B-2a : Comparison Of Experimental Data with 


Joint Failure CDF For Independent Case (15 Deg.) 
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At 8 = 15°, even when the load is closely aligned with the fiber angle, Figure 
B-2a shows that the fiber strength statistics has no effect on the combined 
failure probability by noting that the fiber failure curve is off-scale for the high p 
region. That is, the failure is always dominated by a combination of matrix 
transfer strength and the matrix shear strength. The interaction of transverse 
and shear can be better observed in the linearized Weibull plot, Figure B-2b, 
from which it can be seen that in the lower tail (low p region) the transverse 


Strength dominates and in the upper tail (high p region) the shear strength 








dominates. 
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Figure B-2b : Comparison Of Experimental Data With Joint 
Failure CDF For Independent Case (Linearized, 15 Deg.) 
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While the experimental data are not close to the predicted joint failure CDF 
under the simplification of independence, its significance can not be evaluated 
because of the small sample size. With the eight samples in the current case, 
neither the upper tail nor the lower tail can be observed,thus the failure 
probability coupling remains unresolved. The location discrepancy perhaps 
has more significance. It Suggests that combined transverse and shear 
Stresses have weakened the composites and that the failure mechanistic 


coupling needs to be treated. 
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Figure B-3a : Comparison Of Experimental Data With 
Joint Failure CDF For Independent Case (20 Deg.) 
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At @ = 20°, Figures B-3a and 3b show that even with a slight angle change of 
5°, the effect of the shear component diminishes rapidly. The transition region 
has shifted higher on the upper tail. Physically, it means that unless a large 
number of samples are tested, only the strongest of the samples will fail in 
shear, and the remaining will fail by transverse stress. For this reason, and 
again because of the small number of samples, the shape indicated by the 
experimental data has no significance. The location difference between the 


data and prediction suggests that mechanistic coupling needs to be included. 
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Figure 3-3b : Comparison Of Experimental Data WithJoint 
Failure CDF For Independent Case (Linearized, 20 Deg) 
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At 8 = 30°, Figures B-4a and 4b show that the effect of shear is shifted to 


even a higher upper tail. All comments on 6 = 20° data apply to 8 = 30° as well. 
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Figure B-4a : Comparison Of Experimental Data With 
Joint Failure CDF For Independent Case (30 Deg.) 
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Figure B-4b : Comparison Of Experimental Data With Joint 
Failure CDF For Independent Case (Linearized, 30 Deg.) 
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At @ = 45°,Figures B-5a and 5b show that the effect of shear is practically 


nonexistent; in fact, it can not be observed from the range of the scale 


presented. What is significant is the location of data is below the location of the 


joint CDF suggesting mechanistic coupling. 
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Figure B-5a : Comparison Of Experimental Data With Joint 
Failure CDF For Independent Case (Linearized, 45 Deg.) 
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Figure B-5b : Comparison Of Experimental Data With Joint 
Failure CDF For Independent Case (Linearized, 45 Deg.) 
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At 8 = 60°,Figures B-5a and 5b show that the effect of shear is no longer 


present and that the fit of the data is much improved. The latter observation 


further substantiates the existence of mechanistic coupling. 
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Figure B-6a : Comparison Of Experimental Data With 
Joint Failure CDF For Independent Case (60 Deg.) 
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Figure B-6b : Comparison Of Experimental Data With Joint 
Failure CDF For Independent Case (Linearized, 60 Deg.) 
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From the experimental data available in the open literature, it was observed 
that a much larger number of samples are required to identify the probabilistic 
coupling and dependency of strength. What is evident is that mechanistic 
coupling needs to be included as an extension to the investigation herein. 
Furthermore, experimental design using the probability failure criterion is 


mandatory to optimize experimental testing. 
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